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Abstract 

A dark energy model with a phantom scalar field, an usual scalar field and the string 
field theory inspired polynomial potential has been constructed. A two-parameter set of 
exact solutions to the Friedmann equations has been found. We have constructed such 
stringy inspired potential that some exact solutions correspond to the state parameter 
w DE > — 1 at large time, whereas other ones correspond to w DE < — 1 at large time. We 
demonstrate that the superpotential method is very effective to seek new exact solutions. 
We also present a two-fields model with a polynomial potential and the state parameter, 
which crosses the cosmological barrier infinitely often. 

1 Introduction 

One of the most important recent results of the observational cosmology is the conclusion that 
the Universe is speeding up rather than slowing down. The combined analysis of the type la 
supernovae, galaxy clusters measurements and WMAP (Wilkinson Microwave Anisotropy Probe) 
data gives strong evidence for the accelerated cosmic expansion [U E] • 

The cosmological acceleration suggests that the present day Universe is dominated by a 
smoothly distributed slowly varying cosmic fluid with negative pressure, the so-called dark en- 
ergy (DE) 0, HI [5j El El [8] . To specify a component of a cosmic fluid one usually uses a phe- 
nomenological relation between the pressure p and the energy density g corresponding to each 
component of fluid p = wg. The function w is called as the state parameter. Contemporary 
experiments [U [21 [31 0] give strong support that the Universe is approximately spatially flat and 
the DE state parameter w DE is currently close to —1: 

w DE = - 1 ± 0.2. (1) 

The state parameter w DE = — 1 corresponds to the cosmological constant. From the theoret- 
ical point of view (see [8, 9J and references therein) this domain of w DE covers three essentially 
different cases: w DE > —1, w DE = — 1 and w DE < —1. From the observations there is no barrier 
between these three possibilities. Moreover it has been shown in (TO] that the state parameter 
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w dei which gives the best fit to the experimental data, evolves from w DE ~ to w DE ^ — 1 
and for a large region in parameter space an evolving state parameter w DE is favoured over 

W DE = - 1 - 

The standard way to obtain an evolving state parameter w DE is to include scalar fields into 
a cosmological model. Under general assumptions within single scalar field four- dimensional 
models one can realize only one of the following possibilities: w DE ^ — 1 (quintessence models) 
or w DE ^ — 1 (phantom models) [TT]. Two-fields models with the crossing of the cosmological 
constant barrier w DE = — 1 are known as quintom models and include one phantom scalar 
field and one usual scalar field. Note that the most of phenomenological models describing the 
crossing of the cosmological constant barrier [12j [131 E] use a f ew scalar fields or a modified 
gravity. 

Nowadays string and D-brane theories have found cosmological applications related to the 
acceleration of the Universe. In phenomenological models, describing the case wbe < — 1> all 
standard energy conditions are violated and there are problems with stability at classical and 
quantum levels (see [51 [151 EH] and references therein). Possible way to evade the instability 
problem for models with wde < — 1 is to yield a phantom model as an effective one, which arises 
from more fundamental theory with a normal sign of a kinetic term. In particular, if we consider 
a model with higher derivatives such as (j)e~ n (j), then in the first nontrivial approximation we 
obtain </>e~ D </> ~ <p 2 — </>□</>, and such a model gives a kinetic term with a ghost sign. It turns 
out, that such a possibility does appear in the string field theory (SFT) framework [17] (see 
also [TH [T6]). namely in the theory of fermionic NSR string with GSO— sector. According 
to Sen's conjecture (see [18] for review), the scalar field <fi is an open string theory tachyon, 
describing the brane decay. The four dimensional gravitational model with a phantom scalar 
field is considered as a string theory approximation, that gives a possibility to solve instability 
problems. 

In this paper we consider a SFT inspired gravitational model with two scalar fields and a 
polynomial potential, which is a generalization of a one-field cosmological model, describing 
in [9]. The first two- fields generalization of this one-field model has been proposed in [19] as a 
polynomial model, which has a one-parameter set of exact solutions with the state parameter 
wde, which crosses the barrier w DE = —1 at large time and reaches —1 from below at infinity. 
In this paper we construct a new model with a two-parameter set of exact solutions, for some 
values of parameters we obtain w DE < — 1 at large time, whereas for other w DE > — 1 at large 
time. Note that the different behavior of w DE at large time corresponds to one and the same 
potential and asymptotic conditions of the fields. 

We study different possibility to use the superpotential method and demonstrate that it is 
very useful not only to construct potential for the given exact solutions, but also to seek new 
exact solutions. To demonstrate that the superpotential method allows to find a form of a 
polynomial potential and solutions for the given Hubble parameter we construct a toy two-fields 
model for the Hubble parameter proposed in the SFT inspired model with high derivatives ]14j. 

2 String Field Theory Inspired Two-Fields Model 

We consider a model of Einstein gravity interacting with a single phantom scalar field <fi and 
one standard scalar field £ in the spatially flat Friedmann Universe. In typical cases a phantom 
scalar field represents the open string tachyon, whereas the usual scalar field corresponds to the 
closed string tachyon [T71 [191 ISHl IZE] . Since the origin of the scalar fields is connected with the 
string field theory the action contains the typical string mass M s and a dimensionless open string 



2 



Mo „ 1/1 



coupling constant g a : 

S= Jd'x^ (j^R +^^(8^-8^0 -V(M)), (2) 

where Mp is the Planck mass. The Friedmann metric g^ u is a spatially flat: 

ds 2 = — dt 2 + a 2 (t) {dx\ + dx\ + dxfj , 

where a(t) is a scale factor. The coordinates (t,Xi) and fields and £ are dimensionless. 
If the scalar fields depend only on time, then the equations of motion are as follows 

9y 

+ 3/70=—, (5) 

ay 

^ + 3^=-^. (6) 

For short hereafter we use the dimensionless parameter m p : m 2 = g 2 Mp/M 2 . Dot denotes the 
time derivative. The Hubble parameter H = a(t)/a(t). Note that only three of four differential 
equations ([S])-© are independent. Equation ([6]) is a consequence of ©-([5]). 
The DE state parameter can be expressed in terms of the Hubble parameter: 

2 H 

W DE = "I " ~J^. (7) 

The crossing of the cosmo logical constant barrier Wde = —1 corresponds to change of sign of 
H. The phantom like behavior corresponds to an increasing Hubble parameter. If we know the 
explicit form of fields 0(t) and and do not know the potential V(0, £), then, using eq. PJ, 
we can obtain H(t) with an accuracy to a constant: 

/ t t 

Hit) = ^\ fi 2 {r)dr- U 2 {r)dr\+C. 

At the same time if we know H(t) we can find the potential as a function of time: 

V(t) = ml (3H(t) 2 + H(t)) . (9) 

The Aref'eva DE model [T7| (see also P, El dSl E21 E3] ) assumes that our Universe is a slowly 
decaying D3-brane and its dynamics is described by the open string tachyon mode. To describe 
the open string tachyon dynamics a level truncated open string field theory is used. The notable 
feature of such tachyon dynamics is a non-local polynomial interaction [181 12H [25], [26j [2Tj, 28j. 
It has been found that the open string tachyon behavior is effectively modelled by a scalar 
field with a negative kinetic term [25]. In this paper we consider local models with effective 
potentials V(0, £). The form of these potentials are assumed to be given from the string field 
theory within the level truncation scheme. Usually for a finite order truncation the potential 
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is a polynomial and its particular form depends on the string type. The level truncated cubic 
open string field theory fixes the form of the interaction of local fields to be a cubic polynomial 
with non-local form-factors. Integrating out low lying auxiliary fields one gets the fourth degree 
polynomial [24, 25J. Higher order auxiliary fields may change the coefficients of lower degree 
terms and produce higher degree monomials. 

The back reaction of this brane is incorporated in the dynamics of the closed string tachyon. 
The scalar field £ comes from the closed string sector, similar to [30J and its effective local 
description is given by an ordinary kinetic term [21] and, generally speaking, a non-polynomial 
self-interaction [3 lj . An exact form of the open-closed tachyon interaction is not known and, 
following |19j . we consider the simplest polynomial interaction. 

More exactly we impose the following restrictions on the potential V(<f>, £): 

• the potential is the sixth degree polynomial: 

6 6-A; 

mo = EE c ^> ( 10 ) 

fc=0 j=0 

• coefficient in front of the fifth and sixth powers are of order and the limit — > oo 
gives a nontrivial fourth degree potential, 

• the potential is even: V(4>, £) = V(— <fi, — £). It means that if k + j is odd, then Ckj = 0. 

From the SFT we can also assume asymptotic conditions for solutions. To specify the asymp- 
totic conditions for scalar fields let us recall that we have in mind the following picture. We 
assume that the phantom field <p(t) smoothly rolls from the unstable perturbative vacuum (<fi = 0) 
to a nonperturbative one, for example = 1, and stops there. The field £(£) corresponds to 
close string and is expected to go asymptotically to zero in the infinite future. In other words 
we seek such a function <p(t) that 0(0) = and it has a non-zero asymptotic at t — > +oo: 
0(+oo) = A. The function should have zero asymptotic at t — > +oo. At the same time we 
can not calculate the explicit form of solutions in the string field theory framework. 

In this paper we show how using the superpotential method we can construct a potential 
and exact solutions, which satisfy conditions obtaining in the SFT framework. 

3 The Method of Superpotential 

The gravitational models with one or a few scalar fields play an important role in cosmology and 
theories with extra dimensions. One of the main problems in the investigation of such models 
is to construct exact solutions for the equations of motion. System with a polynomial 

potential V(<j), £) is not integrable. Moreover we can not integrate even models with one scalar 
field and a polynomial potential. 

The superpotential method has been proposed for construction of a potential, which corre- 
sponds to the exact solutions to five- dimensional gravitational models [32] ■ The main ideas of 
this method are to consider the function H(t) (the Hubble parameter in cosmology) as a function 
(superpotential) of scalar fields and to construct the potential for the special solutions, given in 
the explicit form. Let 

H(t) = W(<f>(tU(t)). (11) 
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Equation (j3J) can be rewritten as follows 



d(j) <9£ 2m 2 v 
If one find such W(4>,£) that the relations 



dW ; dW ■ 1 /- 2 



2m ^' ( 13 ) 



2 <^, (14) 




I/ = 3^ + 2^[(^-) -(^) ] (15) 

are satisfied, then the corresponding 4>(t), £(t) and H(t) are a solution of system ©-(E]). 

The superpotential method separates system (jS])-® into two parts: system (fT3l) - (|T4l) . which 
is as a rule integrable for the given polynomial W ((/>,£) and equation (TT5l) . which is not inte- 
grable if V^(0,£) is a polynomial, but has a special polynomial solutions. The way to use of 
superpotential method does not include solving of eq. (ITS]) . The potential V(0,£) is constructed 
by means of the given W(<fi,t;). 

There are a few ways to use the superpotential method. The standard way [32] is to construct 
the potential for the solutions given in the explicit form. We assume an explicit form of solutions, 
find the superpotential W and use ( JTBT) to obtain the corresponding potential V. If we consider 
one-field models, putting, for example, £ = in (J3j) — (JH1) , then from ( |T3|) we obtain W(<f>) up to a 
constant. At the same time solving ( {TBI) we obtain a one-parameter set of solutions: (j)(t — to)- 
So in the case of one-field models we have the following correspondence 

<f>{t-t ) <-> W{<f>) + C, (16) 

where to and C are arbitrary constants. 

In two-fields models the correspondence (fT6|) does not exist and the superpotential method 
gives a possibility to find new solutions. Indeed, equations (|T3|) - (|T^|) form the second order 
system of differential equations. If this system is integrable then we obtain two-parameter 
set of solutions. To assume some explicit form of solutions means to assign a one-parameter 
set of solutions. The superpotential method allows to generalize this set of solutions up to 
two-parameter set. On the other hand we can construct different forms of superpotential and 
potential, which correspond to one and the same one-parameter set of solutions. 

The idea to consider the Hubble parameter as a function of scalar fields and to transform 
([3])-([6]) into ffT3l - (fT5l) has been used in the Hamilton- Jacobi formulation of the Friedmann 
equations [231 E3] (see also [35]) and does not connect with supersymmetric and supergravity 
theories. At the same time the method, based on the idea to apply system (|13|) - (|15p instead of 
the original equations of motion for the search exact special solutions, is actively used in two- 
dimensional fields models [36j [37] and supergravity [38]. Equations (|T3l - (jT4l) are known as the 
Bogomol'nyi equations [39j (see also [37] ). The superpotential method is a combination and a 
natural extension of these two methods. This method is actively used in cosmology [51 \W \ HI] . 
Let us note generalizations of this method on the equations of motion, describing the close 
and open Friedmann universes [JO], systems with the cold dark matter [JT] and the Brans- 
Dicke theory [32] ■ The idea to consider the Hubble parameter as a function of scalar fields 
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and to transform ([HD-© into (|T3l) - (|T5l) has been used in the Hamilton- Jacobi formulation of 
the Friedmann equations [32J El] (see also [22]) and does not connect with supersymmetric 
and supergravity theories. At the same time the idea to apply system (TT^- lfloT) instead of the 
original equations of motion and to seek in such a way exact special solutions is actively used in 
two-dimensional fields models [36, 37J, supergravity [38J and supersymmetric models with the 
BPS states. Equations ([TBI and ffT4"|) are known as the Bogomol'nyi equations [39] (see also [37]). 
The superpotential method is a combination and a natural extension of these two methods. At 
present the superpotential method is actively used in cosmology P d2l EH HI]- Let us note 
generalizations of this method on the equations of motion for the close and open Friedmann 
universes |40j . systems with the cold dark matter [JT] and the Brans-Dicke theory [42J. 



4 The construction of potentials for the given solutions 
4.1 Non-polynomial potential 

In this section we demonstrate that one and the same solutions can correspond to both polyno- 
mial and non-polynomial potentials. In the next section we show that the superpotential method 
allows to find different exact solutions, which correspond to the different behavior of the Hubble 
parameter, but one and the same potential. 

From the asymptotic conditions we assume the following explicit form of solution: 



(Pit) = Atanh(cut) and f(t) = A V 2 0- + ty (17) 

cosh(u;t) 



where A > 0, u > and b > — 1. 
From (1H1) we obtain 



Hit) = — % (3 tanhiut) - (3 + 2b) tanh 3 ^)) . 
6rrip V / 



Note that this kink-lump solution is a natural generalization of the kink solution for the 
one-field phantom model [9]. The behavior of the Hubble parameter at large time depends on 
the parameter b. From the contemporary experimental data it follows that the present date 
Universe is expanding one that corresponds to H > at large time. The condition lim H > 

t=+oo 

is equivalent to b < 0. Eventually, we state that — 1 < b < 0. 

On the other hand, in the past there were eras of the accelerated and decelerated expanding 
Universe, it means that the Hubble parameter H has to be not a monotonic function and should 
has an extremum at some point t c > 0. From ffl8l) we obtain that 

1 , / J2(b + l)(2b + 3)' . 

t c =- arccosh ± v v n J \ J . 19 
uj \ 2(6+1) 1 

We have assumed that b > —1, so the sign "+" corresponds to real t c . At t > the Hubble 
parameter H has one extremum, namely a maximum. The corresponding DE state parameter 
w DE is given by 

2 (2(6 + l)cosh(^) 2 -3-26)cosh(^) 2 
w DE - -1 + 12m p y4 2 sinhM)2(26sinhM 2_3 )2 • 
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It is easy to see that at large time w^e > — 1, so we obtain the quintessence like behavior of the 
Universe^. 

Let us construct a potential, which corresponds to fields The functions <j)(t) and £(t) 

are solutions of the following system of differential equations: 




e 



(22) 



2(1 + b) A 2 ' 

The straightforward use of the superpotential method gives 



W = JlL(A-±dA 9W - ^ 1 g (23) 

d<f) 2m 2 p \ A V /' <9£ 2m2y 2(1 + 6)A 2 ' 1 ; 

Therefore, 



* = W = " T " V 2(1 + 6)^ + ^ I ' (24) 

where Hq is an arbitrary constant. Different values of Hq correspond to different V(4>,£). The 
obtained potentials 

K _^(^_| + _il_ |+3mX (25) 

are polynomial ones only in the flat space-time {m 2 = oo) and do not satisfy conditions of 
Section 2. 

The goal of this paper is to construct a polynomial potential model with such set of exact 
solutions that the quintessence large time behavior corresponds to some solutions and the phan- 
tom large time behavior corresponds to other ones. The potential and solutions should satisfy 
conditions from Section 2. In other words our model should be the SFT inspired one. We make 
this construction in two steps. At the first step we construct polynomial potential for ([17]) . At 
the second step we find new solutions for the obtained polynomial potential. 

4.2 New polynomial potentials for the given solutions 



Let us construct for the functions ffTTj) such a superpotential that the corresponding potential 
has the polynomial form. Functions ffT7|) satisfy not only system (|22l) . but also the following 



x It has been shown [23] that if we consider the other pair of scalar fields 



0(t)=tanh(t), i(t)= ^tff , (21) 

cosh(2t) 

then for some values of parameter b, for example b = —0.01, the corresponding Hubble parameter has both a 
maximum and a minimum at t > and increases at large time. Note that the polynomial potential, which 
corresponds to solutions (|2ip . is not known. 
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system of differential equations: 



Aub 



A 2 



2 - 1 ) + < 



2A' 



The corresponding Hubble parameter (superpotential) is given by 



H = W 



0J(j) 

2m 2 



Ah 



3A 2 



1 + 



2A 



+ H . 



To obtain even potential we put H = 0: 



cj 2 , 2 . 2 3tu 2 2 , „ 
2^ + l^" ' 



3A 2 



e 2 
+ — 

2y4 



(26) 



(27) 



(2? 



This example shows that the same functions <fi(t) and £(£) can correspond to essentially 
different potentials V(<f>,£). So, we conclude that in two- fields models one has more freedom 
to choose the potential, without changing solutions than in one-field models. Moreover, the 
solutions do not change if we add to the potential V (or V) a function SV, which is such that 
SV, d(SV)/d<p and d(SV)/d£ are zero on the solution. For example, we can add 



6V = K(<j>,0 



+ 



2(1 + 6) 



e-A 2 



(29) 



where K(£, <ft) is a smooth function. So, we can obtain new potentials, which correspond to the 
given exact solutions (JT7j) . without constructing of new superpotentials. 



5 Construction of new solutions via the superpotential 
method 

In previous section we have shown how we can choose potential for the given solutions. In this 
section we demonstrate the possibility to find new exact solutions (may be in quadratures) using 
superpotential method. Let us consider the model with the potential (|28l) . It is easy to see that 
system ( 1261) has not only solutions (TPT1) . but also the trivial solutions {</>(£) = iA, £(£) = 0} 
and solution 

(f>(t) = -Atanh(tub(t - to)), f (t) = 0. (30) 

If ^ 0, then, using the second equation of (1261) . we obtain the second order differential 
equation in £(t): 

m= s m) .^ + !lzgn. (31) 

The solutions of eq. (l3Tj) with b > — 1 are defined in quadratures 

t- t0 =±r , =rf6 (32) 

J u^/2A 2 i 2b + 2A 2 b£ 2b + £ 2b + 2 + 2A 2 C + 2A 2 bC 
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where C and to are arbitrary constants. For some values of parameter b the general solution to 
(I26p can be written in the explicit form, for example at b = — 1/2 we obtain: 



<j>(t) 



A ((C?Cf + AA 2 ) - Cfe-^) 



(CfCi + AA 2 ) + 2C\C 2 + C\e-^ 

Ad A 

(C 2 C 2 + AA 2 ) e"" + 2C 2 C 2 + C 2 e~ 



(33) 



■Ult 



where C\ and C 2 are arbitrary parameters. It is easy to check that for all values of Ci, but 
C\ = 0, and C 2 solutions (1331) and the Hubble parameter H(t) satisfy the following asymptotic 
conditions: 

A 2 u 

4>(±oo) = ±A, f (±oo) = 0, H(±oo) = ± — -. (34) 

So we have constructed a gravitational model with a two-parameter set of exact solutions. 
The potential and solutions satisfy conditions, imposed by means of the string field theory (see 
Section 2). 

Let us analyze the property of the obtained solutions and the cosmological consequences. 
System fl26|) is invariant to change £(t) on —£(£), so each solution (f>(t) corresponds to two 
solutions ±£(t). Note that the function <ft(t) is invariant to the change C\ — > —C\, whereas the 
function £(t) changes a sign. The Hubble parameter depends on £ 2 , so, without loss of generality, 
we can put G\ > 0. 

System J26l) is autonomous one, so if there exists a solution {4>{t) , ^(t)} , then a pair of 
functions {(p(t — to),£,(t — t )}, where to G C, has to be a solution as well. It is convenient to 
use in ( )33|) such parameters that one of them corresponds to a shift of solutions in time. We 
put C\ = exp(t ). Using the restriction to G M we come to the condition C\ > 0. For short we 
introduce the new parameter C = C\C 2 instead of C 2 . Solutions (1331) take the following form: 

_ A ((C 2 + 4A 2 )e^-'°) - e -"('-'o)) 



(C 2 + 4A 2 )e"('-'°) + 2C + e-^-^ 
AA 



(35) 



(C 2 + 4A 2 )e<"(*-*<>) + 2C + e - w (*-*o) ' 



To compare the obtained solutions with the initial solution ffl7|) . we introduce new parameter 
t 1 = t + too, where 

t 00 = -^\n(C 2 + AA 2 ). (36) 
Now functions 4>(t) and £(£) are 



<P(t) 



C ^ v / C 2 +4A 2 T C 



(37) 



Let us consider solutions with ti = 0. It is easy to see that in this case 



Aujy/C 2 + AA 2 i / \ . . 

= 0, = ^^=^= > and f (0) = 0. 38 

K J K ! C + VC 2 + AA 2 K ! K J 
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From (j4J) it follows that H(0) > and from (I27p it follows that H(0) = 0. Therefore, solutions 
with t\ = and an arbitrary C are cosmological bounce solution (see, for example, [43]), in 
other words, ait) has a bounce in the point t = 0. 

Let us consider how the behavior of the Hubble parameter H(t) depends on C. 

In the case C = we have solutions 



<p (t) = Atanh(u;(t - £i)) and f (*) 



A 



cosh(u;(i — ti)) 

At £i = these solutions coincide with solutions (TTTI) . The corresponding Hubble parameter 

A 2 cu 



(39) 



# = _^ ^3 tanh(wt) - 2 tanh 3 (cut) ) 



(40) 



has a maximum at the point t max = — In (v2 — l) /a; ~ 0.881/co> and the quintessence large time 
behaviour. The solutions </>o and £q, the Hubble parameter H and the state parameter w DE are 
presented on Figured] (we put A = 1, to = 1 and m 2 = 1/6). 




2 

# 

1,5 



0,5 



1 2 3 4 5 6 




1 2 3 4 5 6 

L 



DE 



Figure 1: The fields <ft an d £ (left), the Hubble parameter if (center) and the state parameter 
wde (right) at C = and £i = 0. 



For arbitrary C the Hubble parameter is as follows: 

A 2 uj (e w (*°-*) + (C 2 + 4A 2 ) e^*-* )) 



6m 2 ( eW (to-t) + 2C + (C 2 + AA 2 ) 



Y e 2w(*o-t) _|_ 6(7e a;< - to ~ t ' ) + 



(41) 



+ 10 (C 2 + AA 2 ) + (C 2 + AA 2 ) e^-^ + (C 2 + AA 2 ) 2 e 2 "^ 

The straightforward calculations give that for all C, but C = ±2A, H{t) = at four points 

I , / AA 2 + C 2 ± 2A\/%A 2 + 2C 2 



t n In 



(C±2A)(C 2 + 4A 2 ) 



fc = l,...,4, 



(42) 



where two signs "±" are independent. Note that if C ^ ±2 A, then H(t m ) ^ 0. Therefore, the 
Hubble parameter H(t) has extrema at points t m . 

At C > 2A all four points t m ^ do not belong to real axis. 

If C — 2A, then H(t) = at two points, which do not belong to real axis: 



t = t - - ln(-2A) and t =t -~ ln(-4A). 

1 oo '■■> 



to 



(43) 



10 



So, at C ^ 2A the Hubble parameter H(t) is a monotonically increasing function and its 
behavior is close to the behavior of the Hubble parameter in one- field model [9]. 

At —2A < C < 2A the function H(t) has extrema at two points. If t\ = 0, the is 
an odd function, whereas £(t) is an even one. Therefore the corresponding Hubble parameter, 
calculated by means of fl27|) . is an odd function. It is easy to check that on semi-axis t > the 
Hubble parameter H(t) is positive and, hence, has a maximum at C < 2 A (see Figure [2]). Thus, 
the behavior of H(t) in the case — 2A < C < 2A looks like its behavior at C — 0. 





1 2 3 4 5 6 



J DE 



Figure 2: The fields cf) and £ (left), the Hubble parameter H (center) and the state parameter 
Wbe (right) at C — 1 and t 1 = 0. 



If C 



-2A, then H(t) = at two points: 
1 



t m _ = t - - \n(2A) 

UJ 



and 



t --ln(4A). 

UJ 



(44) 



At these points H = ±16A 2 /m^ ^ 0, hence, the Hubble parameter behavior is close to H(t) on 
Figures [T] and [2J 

Let us consider the case C < —2A. All four points of extremum (14*2!) are real. It means, that 
at C < —2A we obtain a qualitative new behavior of the Hubble parameter. 

If ti =0, then, as it has been noted above, the Hubble parameter is an odd function. The 
derivative of the Hubble parameter at zero point is positive, hence, H(t) has maximum at some 
point t mi > 0, a minimum at t m2 > t mi and is a monotonically increasing function at t > t m2 . 
Note that w DE < — 1 at t > t m2 . Thus we have found the exact solutions, which correspond to 
the nonmonotonic function H(t) with phantom large time behaviour (see Figure [3]). Using the 
superpotential method we have obtained that the model with the potential 



V 



OJ 




+ey 




(45) 



has two-parameter sets of exact solutions. Note that the obtained solutions have one and the 
same asymptotic conditions, whereas the behaviour of the state parameter w DE turn out dif- 
ferent. So, we can conclude that at large time both quintessence {w DE > —1) and phantom 
(w DE < — 1) behavior of w DE are possible to obtain from the SFT inspired effective model with 
one and the same polynomial potential. 
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H 

50 



J DE~- 




5 t 6 



Figure 3: The fields and £ (left), the Hubble parameter if (center) and the state parameter 
wde (right) at C = — 5 and t\ = 0. 

6 Two-fields model with a polynomial potential and wde 
crossing the cosmological barrier infinitely often 

From the Cubic Superstring Field Theory I.Ya. Aref'eva and A.S. Koshelev obtained that the 
late time rolling of non-local tachyon leads to a cosmic acceleration with a periodic crossing 
of the cosmological constant barrier [H]. At large time approximation, when an open string 
tachyon <fi = 1 — 5<p and \S(f)\ <C 1, the following Hubble parameter has been obtained: 



H = H + C H e~ Zrt sin(2z/(t - t )) 



(46) 



where Hq, C h , r, v and to are real constants. 

In [H] the authors consider non-local model and the corresponding Friedmann equations. In 
this paper we construct the two-fields local model with the Hubble parameter (H6j) in the case 
v = r. For simplicity we put C H = l/(2m£) and construct solutions, which do not depend on 



nip. In this case 



Using 



mt 



H = — e~ 2ri (2sin(rt) 2 - (sin(rt) - cos(rt)) 1 



we can define the following explicit form of solutions: 

5~4> = — 2\prt~ Tt sin(rt), £ = V2re~ rt ^sin(rt) — cos(rt) 
It is easy to check that if 

5~(f> = —=e~ rt ( cos(rt) + sin(rt) 

V T V 

then 



Let construct the superpotential: 



dW 

d5(f) 



dW 



e= - 

>/2r5<j> + 2r£ 



v2 _ rt . , ,s 
sm(rr), 



y/2r8(j> - 2r£ 



so 



W 



V2r^64> - r£ 2 ) + H 



(47) 
(48) 

(49) 
(50) 

(51) 
(52) 
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The potential is (we put H = 0) 



V = -r 2 (e + H 2 -2V2^+^(V2^-e) ■ (53) 

Thus we obtain the explicit solutions and the fourth degree polynomial potential, which 
corresponds to the Hubble parameter from the SFT inspired model with high derivatives. 

Note that the standard method to construct models with scalar fields for the given behaviour 
of the Hubble parameter is the method, which uses V(<p,£) as a function of time [13J. If we know 
the Hubble parameter H(t), then, using (Q, we obtain V(t) and after that we can attempt to 
find the functions (j)(t) and and the potential V(</>, £). Such method is very effective if at 
least one of derivatives either or is such a function F(V) that a form of F does not depend 
on <p and £. For example, if 

V(0,O = W)e Q «, (54) 

where a is a constant, then 

and (jHJ) is a linear differential equation in £ 

i + 3H(t)i + a (3H 2 (t) + H(t)) = 0. (56) 

This equation allows to find if the Hubble parameter H(t) is known [13]. The superpotential 
method is not so effective to seek potential in the form (154")) for the given Hubble parameter H(t). 
On the other hand if the required form of the potential is a polynomial, the superpotential method 
is not less effective and maybe even more easy to use than the above-mentioned method. 



7 Conclusions 

In this paper we have investigated the dynamics of two component DE models, with one phantom 
field and one usual field. The main motivation for us is a model of the Universe as a slowly 
decaying D3-brane, whose dynamics is described by a tachyon field |17j . To take into account 
the back reaction of gravity we add a scalar field with an usual kinetic term. 

We construct a cosmological model with the SFT inspired polynomial potential V(4>,£) and 
find two-parameter set of exact solutions. This set can be separated into two subset such that 
one subset corresponds to the quintessence large time behaviour, another subset corresponds 
to the phantom one. Note that both subsets have solutions, which satisfy one and the same 
asymptotic conditions and the additional condition 0(0) = 0. 

We also construct two-fields model with the fourth degree polynomial potential, which cor- 
responds to the Hubble parameter, obtained in the SFT framework [13] . In this model the state 
parameter wde crosses the cosmological constant barrier infinitely often. 

In this paper we actively use the superpotential method and show that there are new ways 
to use this method in the case of two fields. We can not only construct potential for the 
given solutions, but also find new solutions. In particular superpotential method allows to 
generalize a one-parameter set of solutions up to two-parameter set. The superpotential method 
allows to separate the initial system of motion equations ©-([6]) into two parts. One part is the 
equation on superpotential ffl5l) . which in general case is not integrable, but for many polynomial 
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potentials has special solutions. Substituting these solutions into the second part (system (fT3|) - 
fll4p ) we obtain a system of ordinary differential equations, which is usually integrable at least 
in quadratures. Note that the systems of the type (TTSI) - (fT4l) are actively investigated both in 
mechanics and in supersymmetry theories with BPS states. So, the superpotential method allows 
to stand out from the system of the Friedmann equations a subsystem, which can be integrable, 
even in the case of a nonintegrable initial system. On the other hand this method allows to 
make such a fine tuning of parameters of the considering gravitational models, for example, a 
choose of coefficients of the potential, that the explicit solutions exist. 

For solutions with H (0) = we obtain that H (0) > 0. Similar solutions are known as bounce 
ones (see, for example, [43]). Note, that bounce solutions has been obtained in the SFT inspired 
higher-derivative models [HJ US] . 
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